Measuring the cosmological bulk flow using the pecuhar velocities of supernovae 
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We study large-scale coherent motion in our universe using the existing Type lA supernovae data. 
If the recently observed bulk flow is real, then some imprint must be left on supernovae motion. 
We perform a Bayesian Monte Carlo Markov Chain analysis in various redshift bins and find a 
sharp contrast between the z < 0.05 and z > 0.05 data. The z < 0.05 data are consistent with 
the bulk flow in the direction (l,b) = (290^31°, 20^32°) with a magnitude of Ubuik = I88I103 km/s 
at 68% confidence. The significance of detection (compared to the null hypothesis) is 95%. In 
contrast, z > 0.05 data (which contains 425 of the 557 supernovae in the Union2 data set) show 
no evidence for the bulk flow. While the direction of the bulk flow agrees very well with previous 
studies, the magnitude is signiflcantly smaller. For example, the Kashlinsky, et al. 's original bulk 
flow result of Ubuik > 600 km/s is inconsistent with our analysis at greater than 99.7% confldence 
level. Furthermore, our best-fit bulk flow velocity is consistent with the expectation for the ACDM 
model, which lies inside the 68% confidence limit. 



I. INTRODUCTION 

One of the most recent puzzles raised by the cosmo- 
logical observational data is the so-called "bulk flow" 
- a coherent motion of a large part of our visible uni- 
verse. Originally, Kashlinsky et al. claimed the de- 
tection of bulk flow using the Sunyaev-Zel'dovich effect 
in the Cosmic Microwave Background (CMB) radiation. 
The bulk flow velocity magnitude between 600 km/s and 
1000 km/s directed toward the point given in galactic 
coordinates by {I, b) = (283 ± 14°, 12± 14°) at scales sur- 
passing 800 Mpc was measured. This result was later 
reinforced in [31 H] by using a compilation of peculiar ve- 
locity redshift surveys. Their analysis showed a bulk flow 
of 407±81 km/s towards (/,6) = (287±9°, 8±6°), which 
is still large enough to be at odds with the ACDM model. 
Other work on the topic includes Refs. [51-ITT]. Though 
the directions of the bulk flow are very consistent, the 
magnitudes of the bulk flow velocities are less consistent 
and vary from analysis to analysis. 

In this paper we study the peculiar velocities of Type 
I A supernovae. If the bulk flow is real, then some im- 
print must be left on supernovae motion. Since the 
recession velocities of the high redshift supernovae are 
large and the magnitude of the bulk flow is moder- 
ate, one can expect the largest contribution to come 
from low redshift supernovae. Our results are consis- 
tent with this expectation. We find the bulk flow in the 
direction {l,b) = (290l3^°, 20132°) ^i^^ a magnitude of 
1881 

103 ^'^Z ^ 68% confidence for the low red- 
shift (i.e. z < 0.05) supernovae. 



II. ANALYSIS 

The effects of Type lA supernovae peculiar velocities 
have been studied for a while [El [13] (for some earlier 
studies see [21 [TS]) because of their importance in cos- 



mology as standard candles. In particular, the peculiar 
velocities can change the luminosity-distance relation- 
ship. The measured (perturbed) luminosity-distance de- 
pends on the original (unperturbed) luminosity-distance 
and peculiar velocities of the observer and the source. 

z = z + (l + z)h ■ {Vs ~ Vo), (1) 
Dl{z) = il + 2h-v,~h-Vo)DL{S)- (2) 

Vo and Vs are peculiar velocities of the observer and the 
source (supernova) respectively, n is the unit vector along 
the line of sight, and points from the observer to the su- 
pernova, z and Dl{z) are the redshift and luminosity 
distance from the Type la Supernova measurements. Dl 
and z are the unperturbed redshift and luminosity dis- 
tance. Usually, the unperturbed frame is considered to 
be the CMB frame, which requires setting Vo = 0. One 
may also notice from Eq. ([2| that the physics is not in- 
variant under the exchange of Vo and Vs . This in principle 
gives the possibility to find out the locally preferred refer- 
ence frame that moves with Vo with respect to the CMB 
frame. 

On large scales the universe can be described by an 
isotropic and homogeneous FLRW metric 

ds^ = -df + a{tf [dr^ + r2(d6l2 + sin^ ddcf^)] . (3) 

Here a{t) is the scale factor, and we assume a flat met- 
ric. The supernovae luminosity distance relation in this 
metric can be written as 



Dl{z) = {l + z)- 



dz' 



Ho Jo ^nMil + z')^ + nA 



(4) 



where c is the speed of light. The luminosity distance 
modulus is 



7715(2) = Slog 
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Dl{z) 
IMpc 
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(5) 
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This is the relationship between the luminosity dis- 
tance modulus, mB{z), and the redshift, z, that we will 
use in our analysis. This relationship depends on two ve- 
locity parameters Vo and Vs ■ One of the ways to include 
the supernovae peculiar velocities is to have a different t/, 
for each individual supernova. However, we adopt here 
a more straightforward approach. We assume that for 
a large number of data points Vg can be treated as an 
average velocity of the whole system, in other words Vg 
is the bulk flow. 

If the local reference frame coincides with the CMB 
frame (after we remove the CMB dipole), then Vo = 0. 
Then the only velocity parameter is t/,, which is the 
bulk flow velocity. The bulk flow velocity as a vector is 
completely described by its magnitude and direction (in 
terms of galactic latitude and longitude). This implies 
that one would need to perform a three-parameter fit to 
the supernova data in order to find a possible bulk fiow. 
In a homogeneous universe, the cosmological parameters 
r^M and are the same at all scales. The peculiar 
velocity may of course be anisotropic, but the ratio be- 
tween the energy density components is not. However, 
it is still instructive to vary the values of the cosmolog- 
ical parameters to check for any possible degeneracies. 
In principle, the uncertainties in the background param- 
eters could result in effectively larger error bars on the 
luminosity distance, degrading the significance of the fit. 
We therefore perform a five-parameter fit in addition to 
the three-parameter fit, varying the values of cosmolog- 
ical parameters. To obtain the level of significance of 
the Vg ^ case versus the null hypothesis Vg — we 
calculate the one-parameter posterior probability on Vg, 
marginalized over all the other parameters in the fit. 



III. RESULTS 

From the earlier work [Hi [15] , one can infer that the 
peculiar flow plays an important role for z < 0.1. Ac- 
cording to [U 12] the bulk flow can extend to even higher 
redshifts up to(2; < 0.3). Since the supernova data are 
readily available for these redshifts, we therefore expect 
the bulk flow to leak into the supernova measurements. 

We analyze the Union2 data set from the Supernova 
Cosmology Project (TB] using the Cosmomc [T7] code for 
implementing a Monte Carlo Markov chain search of the 
parameter space. We use a modified version of the like- 
lihood code released with the Union2 data set, adding 
bulk flow parameters to the fit. We include full system- 
atic errors for the data set, as described in I16| . 

To perform the three-parameter flt to the data we 
flx the values of the cosmological parameters consistent 
with a best-flt from the WMAP 7-year data set [27, 28J, 
assuming flMh^ — 0.1334, w = —1, and a flat uni- 
verse. In both cases, we take Ui,h^ = 0.0226, and make 
use of the HST measurement of the Hubble parameter 
Ho = lOOh km s"^Mpc"^ HHj by multiplying the likeli- 
hood by a Gaussian likelihood function centered around 
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FIG. 1: Bulk flow direction in galactic longitude / and galac- 
tic latitude b for a three-parameter fit over lis, b, and I, for 
redshifts z < 0.05. 



Lil^elihood for Buik Flow Velocity (z < 0.05) 
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FIG. 2: Bulk flow magnitude Va for a three-parameter flt over 
Vs, b, and I, for redshifts z < 0.05. 

h — 0.72 and with a standard deviation a — 0.08. We 
adopt flat priors on I and sin(&). We divide the data 
into two sets, doing fits separately for z < 0.05 and for 
z > 0.05. Figures ([I]) and ([2| show the best-fit direction 
and magnitude of the bulk fiow for the three-parameter 
flt for z < 0.05. 

For the purpose of testing parameter degeneracies, we 
repeat the analysis assuming broad priors on Q.Mh'^ and 
w 

0.01 < Om/i^ < 0.99 (6) 
-2 < w < -0.2. 

We note that this is an unduly pessimistic assumption: 
inclusion of other data sets such as Baryon Acoustic Os- 
cillations or CMB would constrain these parameters to 
much narrower regions. However, we are interested in 
investigating the effect of parameter degeneracies, and 
therefore allow the parameters to vary over a wide range. 
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Best-fit Bulk Flow Direction (z < 0.05) 
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FIG. 3: Bulk flow direction in galactic longitude I and galactic 
latitude b for a five-parameter fit over Vs, b, I, QmIi^ and 
w, for redshifts z < 0.05. Crosses indicate the positions of 
the supernovae used in the fit, and the diamonds mark the 
directions of the CMB dipole and the Great Attractor. 



Likelihood for Bulk Flow Velocity (z < 0.05) 




FIG. 4: Bulk flow magnitude Vs for a flve-parameter fit over 
Vs, b, I, Q,Mh^ and w, for redshifts z < 0.05. 



Figures ([s]) and Q show the best-fit direction and mag- 
nitude of the bulk flow for the five-parameter fit with 
z < 0.05. 

We first note, by comparing three- and five-parameter 
fits, that including ^Im and w does not add significant 
degeneracy to the system. This is reasonable since it is 
unhkely that background parameters will be degenerate 
with parameters resulting in anisotropy. However, this 
is still an important check since the uncertainties in the 
background parameters could in principle degrade the 
significance of the fit. 

The striking feature is a sharp contrast between the 
z < 0.05 and z > 0.05 data, shown in Figs. (5j6|. The 
z < 0.05 data are consistent with the bulk flow in the 



FIG. 5: Bulk flow direction in galactic longitude / and galac- 
tic latitude b for a three-parameter flt over lis, b, and I, for 
redshifts z > 0.05. 



Likelifiood for Bulk Flow Velocity (z > 0.05) 
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(29OI31 , 2OI32 ) with the magnitude 
km/s at 68% confidence. In contrast. 



FIG. 6: Bulk flow magnitude Vs for a three-parameter flt over 
Vs, b, and I, for redshifts z > 0.05. 



z > 0.05 data (which contains 425 of the 557 supernovae 
in the Union2 data set) show no evidence for the bulk 
flow. The significance of detection of 7^ (compared 
to the null hypothesis Vs = 0) is 95% for z < 0.05, but 
gets progressively worse as we add higher redshift data. 
The absence of evidence for a bulk flow at higher red- 
shifts may be a consequence of the fact that supernovae 
data is extremely sparse at large redshifts, and should 
not be viewed as the evidence of absence of the bulk flow 
at high redshifts. If our likelihood function in Fig. [6] 
were peaked at zero, we would be able to conclude that 
the data rule out bulk flow of a given velocity at a given 
confldence level. Instead, our likelihood function is al- 
most flat, so the conservative conclusion is that we see 
neither evidence for bulk flow nor its absence. To settle 
this issue we certainly need more high redshift data. 
The direction of the bulk flow for z < 0.05 agrees very 
well with previous studies pj-|t4j. However, the magni- 
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tude is significantly smaller. For example, Kashlinsky, 
et al.'s original bulk flow result of Ubuik > 600 km/s is 
inconsistent with our analysis at greater than a 99.7% 
confidence level. An interesting question is whether this 
much smaller magnitude is consistent with the expecta- 
tion from a ACDM model. Our best-fit bulk flow velocity 
is around 200 km/s. The theoretical estimate of the rms 
amplitude of the peculiar velocity can be obtained from 
the power spectrum of primordial perturbations P{k) as 



i;rms = 47r7j2/(f^)' J P{k)W^{kr)dk, (7) 

where /(il) — SHoftt, while is the velocity window 
function usually taken to be a Fourier transform of the 
Gaussian window in r-space. For the ACDM model, the 
rms amplitude of the peculiar velocity at some distance 
d is [2] 

v,„,^ = 250(100/1"^ Mpc/d). (8) 

A redshift of z = 0.05 corresponds to c? = 150h^^ Mpc, 
so the expected bulk flow velocity on that scale is about 
170 km/s, which is inside our 68% confidence limit. This 
implies that our findings are consistent with expectations 
from the ACDM model. 



IV. CONCLUSIONS 

Recent results on cosmological peculiar velocities have 
garnered considerable attention. The original study of 
galaxy clusters using Sunyaev-Zel'dovich effect [T] and 
most of the follow up work report an unusually large bulk 
flow which is at odds with the ACDM model. There are 
however results questioning the original analysis [20H22] . 



Any large scale coherent motion of matter in our uni- 
verse must also be reflected in the peculiar motion of 
supernovae. In this paper we performed analysis of 557 
supernovae in the Union2 data set using the Cosmomc 
code, including full systematic errors for the data set. 
For low redshift supernovae (i.e. z < 0.05) we do flnd 
a bulk flow in a direction consistent with previous work, 
however, with a much more moderate magnitude (we find 
Vbulk = 1881 

103 km/s at 68% confidence, instead of the 
claimed 600—1000 km/s). Compared to the null hypothe- 
sis, the significance of detection is 95%. We find no signif- 
icant evidence for bulk flow for high redshift supernovae 
(i.e. z > 0.05). This absence of evidence should not be 
understood as the evidence of absence. The current high 
redshift supernova data are too sparse to place a signifi- 
cant constraint on large-scale coherent motion. We also 
find that the expected bulk fiow velocity of the ACDM 
model lies inside the 68% confidence limit of our best-fit 
bulk fiow velocity. Supernova data were also analyzed 
in this context in though using different data 

sets and/or different methods. The numerical results are 
mostly consistent with ours, though not necessarily the 
conclusions. We can summarize our results as follows: 
The z < 0.05 supernova data are completely consistent 
with a ACDM model and inconsistent with Kashlinsky, et 
al. [T] to high significance (greater than 99.7% confidence 
level). The z > 0.05 data are inconclusive. 
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